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Abstract 

We discuss a general formalism for numerically evolving initial data in general rela- 
tivity in which the (complex) Ashtekar connection and the Newman-Penrose scalars 
are taken as the dynamical variables. In the generic case three gauge constraints 
and twelve reality conditions must be solved. The analysis is applied to a Petrov 
type {1111} planar spacetime where we find a spatially constant volume element to 
be an appropriate coordinate gauge choice. 
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1 Introduction 



Since their introduction by Abhay Ashtekar in 1986 Ashtekar variables have been the 

basis of remarkable conceptual and technical advances toward a quantum theory of gravity ji) 
Q . (Briigmann maintains an updated e-mail bibliography on publications related to Ashtekar 
variables ||] . ) Applications have also been found in classical general relativity. The most fruitful 
of these have probably been in the study of Bianchi cosmologies . It was suggested 

soon after the formalism was discovered that it might lend itself to classical numerical work, 
mainly because of the simplicity of the constraints in terms of the new variables. Meanwhile 
Capovilla, Dell, and Jacobson (CDJ) discovered a procedure for solving algebraically all of the 
diffeomorphism constraints . This work has opened up new connection-based approaches to 
general relativity. (For further references and an informative schematic tracing the roots of this 
current industry to the work of Plebanski [^, see the review by Peldan |l2j.) Our intention 
here is to introduce the CDJ formalism in a general form suitable for undertaking a finite 
differencing approach. The program will then be applied in a testbed which has already served 
the numerical relativity community well in developing stable and accurate one-dimensional 
code: planar cosmology. Recently, Miller and Smolin have proposed a discretization scheme 
which is based on the same (CDJ) action jl^. The variables and discretization they choose are 
likely most suitable for quantum lattice calculations. 

Our plan is first to introduce Ashtekar variables, displaying the full constraint, evolution, and 
reality conditions in the 3-1-1 formalism. In the next section we derive the general dynamical, 
constraint, and reality relations using a variation of the (CDJ) action in which a complex S0(3) 
connection and a symmetric complex second rank SO (3) tensor are our variables. These latter 
objects turn out to be linear combinations of the five complex Newman-Penrose scalars, and 
we present arguments why we might wish to specify them freely and to monitor their evolution. 
In Section 3 we apply the formalism to a spacetime which possesses two commuting spacelike 
Killing fields. The simplest such model requires only two real Newman-Penrose scalars, and is 
of type {1111}. It is not possible to treat either flat, plane sandwich wave, or colliding plane 
waves (prior to collision) in this formalism. Our variable and coordinate choices are preserved 
in a spatially constant (but time dependent) volume gauge. 

We conclude with a prognosis of ongoing work in solving numerically the planar model pre- 
sented in Section 3 and a discussion of potential future applications. Appendix A presents a 
translation between planar Ashtekar and York variables employed in work by Anninos, Cen- 
trella, and Matzner. In Appendix B we derive the relation between our scalar variables and the 
Newman-Penrose scalars. 



2 Ashtekar Variables 

The constraints and evolution equations in the Ashtekar formalism can be most easily obtained 
through a variation of the Hilbert action amended with a term which vanishes automatically 
by virtue of the cyclic Bianchi identity . We first rewrite the Hilbert action in terms of 

tetrad fields E"j, using the associated Ricci rotation coeficients j and curvature 7?.^^,^ j 

n^J J = 2d[^ii,{j + 2f][^,r!,]^i. (1) 

The Ricci scalar TZ is then just 

n = E'^.Ef'j'R^p" = E^jE^j-R^p^f^'^ . (2) 
The vanishing object which we add is simply —i times the scalar formed from the dual *TZ^^'^ 



1 



Our conventions are that greek letters (/i, i^, p, ...) range over the four spacetime coordinate 
indices, while upper case latin letters from the middle of the alphabet (/, J, K, ...) range over the 
four internal S0(l,3) indices. The S0(l,3) indices are raised and lowered with the Minkowski 
metric, which we take to be rjjj = diag{—l, 1, 1, 1). Lower case latin indices from the beginning 
of the alphabet (a, b, c, ...) range over the three spatial indices, while lower case latin indices from 
the middle of the alphabet (z, j, k, ...) range over the three internal S0(3) indices. Since we shall 
frequently deal with concrete indices, we find it useful to represent contravariant objects with a 
capital letters, and corresponding covariant objects with lower case letters. Thus, we represent 
the covariant tetrad fields as e^^, with inverse contravariant field E'^j. Where possible we shall 
distinguish between four-dimensional objects and three-dimensional objects by representing the 
former in either latin script or capital greek letters. 

We achieve a suitable 34-1 decomposition by taking E^q to be normal to our to our spacelike 
foliation 

E^„^ iN-\-N-^N''), (4) 
where N and TV" are the lapse and shift, respectively. Then 

E^^ = (5) 

where is the triad field. The canonical fields in the Ashtekar initial value formulation of 
general relativity are the densitized triads and the connections A'^^. Let 

t:^dett\. (6) 

The densitized triad of weight one, T°, and the densitized lapse of weight minus one, N, are 
thus 

f'^,:=tTl (7) 

N r^N. (8) 

We can of course express the metric in terms of T", but Rovelli has shown that the densitiized 
triad itself has a simple geometrical interpretation: corresponds to an S0(3)-Lie-algebra- 
valued area two-form, and the S0(3) norm of this two-form is the area [ pT[ . 

The Ashtekar connection is constructed from the 3-dimensional Ricci rotation coefhcients 
uj^''^ and the extrinsic curvature fCab in the following manner. 



where 

and w'^j is the dual of w^'-* , 



A\:=uj\-tK\, (9) 
K\ := T"/Cah, (10) 



1 yfc jfc 



(Note that since the SO (3) indices are raised and lowered with the Kronecker delta, we can 
raise and lower them at our convenience.) 

The curvature -F'^f,'"' associated with the the Ashtekar connection is 

Fab' = 2%A]'^'+2A[/-A,,'=^- 

= 2d^aA,;' ^2A\^A\^, (12) 
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and the amended Hilbert action S takes the form 



Jd'^x iifliA\ - VaA^) - iN'^f\F^l + \Nf-{f^,F:^). (13) 



In (13) we have introduced the covariant derivative operator defined through 

VaA' := daA' + A\ (14) 

where 

:=f]*o-ie*^'=f]o°'=. (15) 

Following the Palatini method, we now vary T", A^ N, iV", and A' independently in the 
action (|l3|). The resulting scalar diffeomorphism constraint C takes the form 

C := f^,f)Fj = 0, (16) 

while the vector diffeomorphism constraint Ca is 

Ca := Fjf\ = 0. (17) 

We have in addition a gauge constraint C which reflects or freedom to perform local triad 
rotations 

C Vafi = 0. (18) 

Our equations of motion are 

A\, = VaA^ + N'F,l + iNf^Fj , (19) 

fa . ^ gufcya _^ 2VbN^^f"'}) + ie'^'^VbiNf^jff). (20) 

Finally, we must restrict ourselves to a real metric. A polynomial form of these reality 
conditions is achieved by requiring first that the doubly densitized contravariant metric g""^ := 

3(Tf ) = 0. (21) 

The second reality condition stems ultimately from our definition ^ of the Ashtekar connection. 
If our triads were real (as we shall assume below in our planar example), we would require that 
the real part of the Ashtekar connection be precisely the Ricci rotation coefficient derived from 
the triads. It is clear that this condition cannot involve explicitly a gauge choice. (We employ 
a generalized notion of gauge which encompasses not only the traditional internal gauge, but 
also coordinate choices. So we are referring to the freedom in N and N"' as well as A^.) The 
required conditions can be obtained most efficiently in polynomial form by requiring that the 
time derivative of g'^'' be real We find that with the choices A"^ — and N'^ = in the 
equation of motion (pO|), 

1'^' = 2^iyey■feI?e(^^/'l)T^'. (22) 
So our second reality condition is 

me^JMf)f%)f'^^ = 0. (23) 

The approach we have described seems at this stage to offer no practical advantage from a 
numerical perspective. Indeed, we may wonder whether any gain at all may be achieved given 
that we have increased the number of constraints by three, and we have twelve non-trivial 
reality conditions in the generic case. Nor do the evolution equations undergo any significant 
simplifications. Ultimately, of course, it must be possible to choose gauge conditions such that 
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the conventional and the Ashtekar numerical approach are exactly equivalent. One merely 
constructs the real part of the Ashtekar connection from the assumed functional form of the 
real triad fields, as prescribed in the definition (||). Thus after fulfilling all reality conditions 
in this manner, the formalism simply amounts to a change of variables, which may or may 
not be useful from the point of view of generating numerical solutions. This point of view is 
illustrated in Appendix A in which we show the equivalence with a numerical construction of 
planar cosmologies due to Anninos, Centrella, and Matzner. 

Thus a truly new method demands the use of new gauge conditions. It is to this strategy 
we turn next. 



3 A Connection Approach 

The Ashtekar approach was itself stimulated by work in complex general relativity, and it was 
from this perspective that a procedure was discovered by Capovilla, Dell, and Jacobson for 
solving the four diffeomorphism constraints ( ]l6| ) and ( |l7| ) pO| ]. They discovered that if is 
an arbitrary traceless, symmetric SO (3) tensor, then the following densitized triad satisfies the 
diffeomorphism constraints 

f\ = S''^.*-!, (24) 
where is the double dual of the Ashtekar curvature 

B". l^^^'^e^.^nj'. (25) 

In hindsight this result is not a surprise since the diffeomorphism constraints are merely alge- 
braic and polynomial in and . Indeed, in substituting our Ansatz ( p^ into the constraints 
(|l|) and ^ we find 

C = eabce'"'f'\f\B\ = m-'^u, (26) 

Ca = eabcB"f\ = tabcB"B''=■^,^, (27) 

where ^ := det{^ij) and B := det{B'^^). We note in passing that the traceless condition could 
be relaxed if we wished to contemplate degenerate metrics, in which case we would permit 
B — 0. However, henceforth we shall assume that B 0. 

The components of are simply linear combinations of the complex Newman-Penrose 
scalars 'i/i, where i = 0, ...,4. In Appendix B we show that 

*ii = i(-*o + 2«'2-*4), (28) 

*12 = ^(«'o-*4), (29) 
*13 = *i-*3, (30) 

*22 = i(^'o + 2«'2 + *4), (31) 

*23 - -^'^'l-^^'3, (32) 
*33 = -2*2, (33) 



with the inverses, 



1-0 = 2(^*11 +*22- 2^^-12), (34) 
*1 = ^(*13 + **23), (35) 
*2 = ^(*ll + *22), (36) 
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^-3 = ^(-*13 + i^'23), (37) 
*4 = ^(-*11+'J'22 + 2^^-12), (38) 

These scalars possess a readily measureable physical significance pointed out by Szekeres 
p3| . ^'o and ^'4 measure transverse contributions to pure gravitational waves, while and 
\l/3 correspond to longitudinal contributions. ^'2 is a Coulombic component which we do not 
permit to vanish since we require that exist. Another quantity of physical interest is 
the Bel-Robinson superenergy tensor whose positive-definite properties can often be useful in 
monitoring gravitational wave disturbances. Its tensorial expression in terms of the Newman- 
Penrose scalars is given in Appendix B. We view these physical interpretations as motivation 
for directly monitoring the dynamics of the Newman-Penrose scalars. We must note, however, 
that in assuming that is invertible, we are excluding Petrov types {31}, {4} and {-} from 
consideration, as we discuss in the Appendix Q ||25|| . 

We shall now derive the equations of motion for ^fy . They can fact be obtained from the 
following action |Q 

S = ijd^x {{A\ - VaA')B'^^-^~\ -f M»e'''*jfe + P'^^^)■ (39) 

(We have added the Lagrange multipliers and p to enforce the symmetry and tracelessness 
of \I'.) However, we shall simply substitute the Ansatz ( p^ for T'J directly into the equations 
of motion (|l|) and 

The characteristic equation for an arbitrary 3x3 matrix '^ij is 

= ^A^^,, - wu^^, + ~ ^M^.. (40) 

Thus in our case 

= - (41) 



n - B-^^-^^.k'fM - ^S.^^ki'^ik). (42) 



Our Ansatz for takes the form 

1 
2 

Substituting into the equation of motion ( [T9| ) we find 

A\ = VaA' + eabcB'''N- + tNeabce'"'B\'f-]^B\ (43) 

= VaA' + eabcB''N' ~^NB^-%,,^,k'^k^■ (44) 

The equation of motion for '^ij requires some more work and we will not include all of the 
details. For this purpose we require the time derivative of 5° 

B", - le'^'^h: (45) 
= e'''"'VbA\. (46) 

So according to our Ansatz ( p^ ) we have 

-ie^^'^VbiNB^-^k'ifM) + ie''''^NB-^-\k-9ki'fi,Vb'^-]^ 

-B^'^~l^rnn^~„\. (48) 
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F.J'^A^ = 2V,,V.A^ (49) 



In (48) we used 

and the evolution equation (^). Comparing the two expressions for T^^, ( pO| ) and (^8|), we find 
ultimately that 

= 2eH(.*,)/^' + A^^^^a^-y - iNVa-^kjeukB'^rr.'^-li- (50) 
In deriving ( |50| ) we made liberal use of the gauge constraint (|l8|). We note here that since 

VaB^, = 0, (51) 
an equivalent form of the gauge constraint is 

C, = B'^Pa-^-], = 0. (52) 
Notice that is manifestly traceless. It is also symmetric, since 

e'^'^^y = -^NVa^kjB^i^'\^+iNVa^nB^i^']^^ 

= iN^iikjE^'iVa^-]^ + 2^N€^uAa^^l^B'',^-]k = 0. (53) 

The first term in the last line vanishes due to the gauge constraint ( p^ ) , while the second term 
is zero because of the symmetry of . 

4 Application to Planar Cosmology 

We shall now apply the preceeding formalism to a spacetime which is presumed to possess two 
spacelike commuting Killing vectors. In this case it is possible to introduce a spatial coordinate 
= z and time t such that the only nonvanishing densitized triad and connection variables are 
A^2^ ^^21 ^^3> K\, K\, K\, and K% [|^. We attempt a further simplification 
motivated by our experience with conventional variables described in Appendix A. We take the 
diagonal components of A^^ to be pure imaginary, and the off-diagonal components to be real, 
so that our non-vanishing connection components are, in addition to the real A^2 ^-nd A^^ (see 
the definition (||)), 

A\ -iK\, (54) 
A^2 := -iK% (55) 
:= -zK%, (56) 
where the K^s are real. The double dual B'^f defined in is therefore ( / is d/dz) 



A\A\ + K\K\ 



B\=\ -iie^-iA\K\ A^^ + K\K% |. (57) 

Let us take \E'y to be diagonal, so that it has the form 

* = I -b I , (58) 
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with a and b real. Thus we have selected the following nonvanishing Newman-Penrose scalars 

vl-o = lia-b), (59) 

^2 = -^ia + b), (60) 

^-4 = lia-b) = *o. (61) 

This corresponds to a spacetime of Petrov type {1111} We note that we can treat neither 
plane sandwich waves nor colliding waves within this formalism since in both cases there exist 
open regions which are of Petrov type {4} and/or {-}. Our spacetime possesses equal left 
and right moving transverse waves, described by and ^'4, respectively, and an everywhere 
present Coulombic part described by ^'2. As mentioned in the introduction, we can now easily 
monitor the Bcl-Robinson super-energy tensor. In Appendix B, we show for example, that the 
super-energy density for observers in the orthonormal basis E'^j is 

Toom^^ia^ + ab + b^). (62) 

Isenberg, Jackson, and Moncrief have investigated the behavior of this object in a Gowdy T^xR 
planar spacetime |p6|| . 

Then our densitized triads (|4|) are 

/ -a-^B\ -b-^B\ \ 
f 'J = -a-^B\ -b-^B\ . (63) 

V {a + b)-^B\ ) 

We shall require that be real. Referring to B";^ in ([57|), we are led to our first rcalitiy 
conditions 

K^^ - A\K\ = 0, (64) 
K\' + A^iK\ = 0. (65) 
Thus our densitized triads are diagonal 

/ -a-^B\ 

f 'J = -b-^B\ I , (66) 

V {a + b)-'B% 

yielding the metric 

{B\)-^B\B%b{a + b) 



'Jab 



{B\)-^B\B\a{a + b) \. (67) 

{B\)-^B\B\ab 



We shall also assume that the only nonvanishing gauge functions are := iA^, N^, and 
N. Substituting T° from (^6|) into ( p^ ) we find the equations of motion 

k\ = -A^A\ + Nab-\a + b)-^B\B%, (68) 

k\ = A^A\ + Nba-^{a + b)-^B\B%, (69) 

k\^A^'+N{ab)-\a + b)B\B^2, (70) 
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A^^^ A^K\~ N'^B\, (71) 
A^^ = -A^K\ + N^B\. (72) 
The remaining equations of motion are (from (^8|) and (|5^)) 

a = N^a' + NK\B\b-\b + 2a) + NK%B%{a - b){a - b)-\ (73) 

b = N^' + NK\B\a'\a + 2b) + NK\B\{a - b){a ~ b)~^ . (74) 

Let us next impose the reahty condition (p3|). In our case this is simply the condition that 
A^2 ^-i^d A^i are the only non-vanishing Ricci rotation coefRcients formed from our diagonal 
densitized triads (p^): 



?3 
3 



(75) 



It turns out that the gauge constraints Ci — Q are now automatically satisfied. This fact is 
easiest to see referring to (|lj) and the connection definition (^: 

C, = daf"-^ + uo^'f''^ ~ le'^'^K^f"^ = 0. (77) 

The sum of the first two terms vanishes (the densitized triad is covariantly constant). The 
second term is the gauge constraint in a real triad formalism in which the canonical variables 
are and K^^^ |^ . In the present case this constraint vanishes identically since both T'\ and 
K^g^ are diagonal. 

Next we must insure that our vanishing variable components remain zero under time evo- 
lution. With our choices of gauge functions, the only non-trivial condition arising from the 
evolution of ^ ij in (|5^) and in ( |l9| ) is 

-i*i2 = {b-a){A^ - N^K\) - Nb'B\{a + b)-^ + NA\B\b-^{a + 2b) = 0. (78) 



Now we must require that the gauge choices (64) and (|65|) are preserved under time evolution. 
We have 

-iB^2 = (-4^ - N^K\)B\- NA^2a^^b-^{a + b)]3\ + {Nba-^{a + b)-^]3\]3^S = 0, (79) 



= -{A^-N^K\)B\-NA^^a-H-^{a + b)B\-[Nab-^{a + b)-^B\B\^ =0, (80) 

Two of the preceeding three relations (|7^), ([t^), and (|80| ) fix the gauge functions A? — N'^K^ 
and iV, so the third must be a new constraint. We present the results without proof as the 

calculation is straightforward but tedious. The constraint turns out to be the requirement that 
the determinant of the densitized triads be spatially constant: 

o.f.,r.nnr = (™i)'. (81) 

We find that TV satisfies the first order differential equation 

TV'Ty-i - A^2B\h-^{a + 2b){a + b){b - a)-^ + b' B\{b - a)-^ = 0, (82) 
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so that — N^K^^ may be expressed most simply as 

j3 _ ^3^3^ ^ -N'{a + by^B%. (83) 

The constant volume gauge has been examined by Rovelli in the context of canonical quanti- 
zation, although he specializes to volume elements which are time-independent [p^ . 

Finally, the shift is fixed through the requirement that the volume clement remain 
spatially constant under time evolution: 

= f ' = 2 (^N^'f + Nf{K\f\ + K\f\ + K\f\)^ ' , (84) 

so that 

iV^" = - (N{{K\f\ + K\f\ + K\f\))', (85) 
and this equation is readily integrable. 



5 Discussion 

We have developed a general formalism for fixing and evolving vacuum initial data consisting 
of a complex SO (3) connection and the set of Newman-Penrose scalars. Although it is certainly 
true that the constraints, reality conditions, and evolution equations are more complicated than 
those obtained using conventional variables, we believe that spacetimes may well exist for which 
this complication is compensated by our ability to fix arbitrarily and then to monitor the time 
evolution of the physically significant Newman-Penrose scalars. 

We list here the full set of relations obtained in Section 3 for the planar case. The dynamical 
variables are a, 6, K\, K\, and K\, with gauge functions N^, N, and A. All are real, and 

they depend only only the coordinates t and z. All constraints are satisfied identically by our 
dynamical variables, but we have the following set of reality conditions: 

K'\ - A\K% = 0, (86) 

K'\ + A\K\ = 0, (87) 

A' -co' - ^ _ a ( aB\B% , 




.2 2 1 T%T\\ b / bB^^B% , 

The equations of motion are 

d = N^a' + NK\B\b-\b + 2a) + NK\B\(a - b){a - (90) 

h = N%' + NK\B\a-'{a + 2b) + NK%B%{a -b){a- b)-\ (91) 

k\ = -A^A\ + Nab-\a + b)-'B\B\, (92) 

k\ = A^A'^ + Nba-'{a + b)-'B\B%, (93) 

k% = A'^ + N{ab)-\a + b)B\B\. (94) 
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The gauge conditions are 



(6 - a){A^ - N^K%) - Nb'B%{a + h)-^ + NA^^j3\h-^{a + 26) = 0, (95) 

{A^ ^ N^K\)B^^- NA^^a~^h'^{a^b)]3\ + {Nba-^ {a + h)-^ B\B^^' = 0, (96) 
{A^ - N^K\)B\ + NA^^a-H-^(a + b)j3\ + [Nab-^{a + b)-^]3\B^.;^' = 0, (97) 

iV^" = - (^iy((ii:\'f \ + iC^a^^l + K\f\)^ ' . (98) 

We are currently finite-differencing these relations. Our strategy is to solve the reaUty 
conditions for the variables K\, and K\ so that the Newman-Penrose scalars a and b 
may be specified freely. We maintain that the complexity of these relations is in itself not 
compelling reason to reject this approach. The issue that counts is whether it is possible to 
construct a stable and accurate numerical code, and the jury is still out. 

We anticipate that it might turn out to be useful to combine the connection approach 
with the conventional one, possibly in choosing one technique for evolution and the other for 
constraints, or applying both to monitor accuracy. We are also investigating the possibility of 
applying the connection approach to the characteristic initial value problem, where the spinorial 
formalism has proved especially useful in analyzing gravitational waves ||25| . 



6 Acknowledgements 

This work was supported by NSF grant number PHY 9211953. D.C.S., A.A., D. M., L.T., and 
B.T. would like to thank the Center for Relativity at The University of Texas at Austin for its 
hospitality. 



Appendix A 



We shall give here the change of variables between our variables T° and -F^'^ and the York 
variables employed by Anninos, Centrella, and Matzner (ACM) in a series of papers in which 
they develop a numerical code for solving Einstein's equations in a plane-symmetric spacetime 
[|16) |18|. The ACM line element is 



ds^ = -a'^dt^ + (jy^dx^ -j- (jy^h'^dy^ + (j)^{l3dt + dzf , (99) 
corresponding to the covariant tedrad field (the components are functions only of z and t) 



a 








fh 





\ 









(100) 



The Ashtekar connection is 

ia-^{2(i(t)(j)' - 200) 
A\ = I -20-V' 




0-i(20'/i-H0/i') 
ia-^{2(j)(l)'l3h + (lPh'13 - 200/1 - 02/i) 




m"H2#'/3 + 0^/3' -200) 
(101) 
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Thus 

/ ~a-^{2[3(t)(j)' ~ 2(j)^) 

K\=i -a-^{2(l)4)'l3h + 4>^h'l3-2(t)^h-(t)^h) 

V -a-HW/3 + 02/?'-20(/)) 

(102) 

with corresponding equations of motion 

4> = (34)' + ]^a(t>-^K\, (103) 

h = a(l3-\K\-K\) + h'f3. (104) 
The remaining equations of motion are 

K\ = A^^A^ + f]{K\' + A^iK%') - a(l)~^h-\hA^i' + hK\K\ + A^^A^^ + K\k\), (105) 

k\ = -A\A' + p{K^^ - A\K%) - ar'{-A\' + K\k\ + A\A\ + K\k\), (106) 
k\ = -A^' - a^-^h-^hA^^' + hK\K\ - A^^ + K\K\). (107) 
The non-trivial constraints are the scalar constraint ([l^): 

A\A\ + if + h{A^; + K\K%) - A^^ + K\K\ = 0, (108) 

and the vector constraint C3: 

/i(7f\' + A^^K\) + K"^^ ~ A\k\ = 0. (109) 

The gauge functions A'^ := iA"^, a, and are partially determined through the conditions 
that f\ =f\=0 and f\ = f%. We find 

A^ = -(3K%+a'(l)-'^, (110) 

and 

(3' = arHK%-K\). (Ill) 

As noted in Section 2, there are no further conditions, since all of our variables are real, and 
we have constructed the only nonvanishing w^^'s from the triads: 

uj\=A\=2(t)-^(t)'h + h', (112) 

uj\ = A^i = -20- V- (113) 

We now give the change of variables to the ACM formalism, where a York conformal de- 
composition is employed. ACM define a traceless extrinsic curvature 

Aab ICab - -^gablC, (114) 

where /C := /Cq- The conformally transformed Aab is 

Aab := 4^ Aab- (115) 
The independent variables are taken to be 

fi:^A\-A\, (116) 
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in addition to A^^ := A and IC. They are obtained from our variables as follows: 

fi = cb\K\-h-'K\), (117) 

A^^cf>\2K%-K\-h-^K\), (118) 

IC ^ (j)-^{K\ + h-^K'\+ K%), (119) 
Substituting into the constraints (|l^) and (^l]) we find for the scalar constraint 

(,/>'/.)' = ^ i^^2h~'h" + ^JC^cj^^ - i(f + 3A')r') , (120) 

while the vector constraint is 

A' - hi-^h'A + hi-^h'fi - '^(t>'^^' = 0. (121) 

As we pointed out in the text, it is by no means obvious that the system of evolution 
equations, constraints, and gauge conditions ( |103[ ) - (|lll| ) are in any way superior from a 
numerical perspective to those obtained through the York procedure. We have written a fully 
constrained code in which we solve the constraint ( |108| ) for K\ and substitute into the vector 
constraint to obtain a first order ordinary differential equation for K\. The metric functions h 
and (j) then chosen to represent incoming colliding waves. Data is evolved using a second- 
order accurate leapfrog procedure with downstream updating. This is, however, a procedure 
which is also available using conventional variables. Furthermore, it could be argued that the 
York procedure possesses greater physical motivation, since one is free to specify a conformal 
class of metrics. Finally, it is much more difficult to implement periodic boundary conditions 
with our choice of independent variables. On the other hand, the connection approach we 
describe in the text is well motivated physically, and there is no difficulty in implementing 
periodic boundary conditions. 



Appendix B 

We shall show that the components of the symmetric SO (3) tensor "^ij used in the construction 
of the densitized triad T°'^ in (^) are the linear combinations of the Newman-Penrose scalars 
given in (^ - 13^). For this purpose it is most convenient to work with both SU(2) and SL(2,C) 
spinors. Although all of the results we shall give here are known, they are not to be found 
summarized in one location, and we think our precedures may also be more accessible to non- 
experts in spinorial calculus. In the course of this derivation the relation between the tetrad 
and triad 3-1-1 decomposition will also be made apparent. 
We let 

^AA' represent times the 2x2 identity and Pauli malices A' 

- 71^^ - ( ! ; ) ' (124) 
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The indices A and A' range over and 1. For visualization purposes we shall conceive of the 
first index as representing the row and the second as representing the column. These SL(2,C) 
indices may be raised and lowered with the Levi-Civita symbols 

CAB = <^A'B' = - 1 ) ■ 

Note that we are letting a capital latin letter from the middle of the alphabet represent a 
Minkowski index ranging from to 3. 



,1 A 



2 A 



,3 A 





A 

3 




' 








^ 


71 ( 


i 




' 1 


71 ( 





1 



(127) 
(128) 
(129) 



We have the following useful identities: 



/ s^J BA' _ '-^iJxB 



S 4 /I'S 



= \v"Si + i6^i6i^ + -j=<5f<5/.^^'=)a'= (130) 

<^i^^<^3 BC = 2'5y^s + ■^'^'J''^'' ^B- (131) 

With the aid of these matrices and the tetrads given in (|^) and (H) , we can introduce the 
soldering forms S''^^/ : 

EV (132) 

As the name suggests, these objects 'solder' the tangent space to spinor space, and we shall 
now use them to effect our translation from tensorial to spinor expressions. For this purpose 
we construct a complex connection ^"^ 



A^," ■.= -{n^"~^*n^''), (133) 



where 



■.^^e^^'^n.KL, (134) 



and e^-^^^ is the Levi-Civita symbol with e"^^'^ — 1. Since 



ij 



A'' = iA'', (135) 



^„ is said to be self-dual. Self-dual objects have natural spinorial analogues. In particular, 



we assert without proof that the spinorial equivalent of the tensorial curvature J^^^ ^'^ is |3C 



T IJ — S^I s^-J ID ASIA'S' _ s^I s^J A''r, AB n 1<i\ 

~ ^ AA'^ BB' '"^iiu £ — ^ AA'^ B 'V"^ ' U'JOj 

where 

T^,^ " := 2d[,A^^ + 2^[^ (137) 

n^^"^^ is the spinorial curvature constructed from the spinorial connection Xt^^^ which leaves 
the metric cab covariantly constant 

X/^ = ■'^'sV'f^, /./. (138) 
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Comparing yl'^ in ( 133 ) with the definition of the Ashtekar connection A\ in (^), we discover 
that 

A\ - Ia\. (139) 

Thus 



T — 98, A "^^ 9 A, ^'^ A ^+9Ar A ^ 

■^ab — '^"[a'^b] + ^-^[a •^b]k + ^"^[a ■^b]0 



- ALA' 



[a^ b] 



2^ab ' 



(140) 
(141) 



where in (|140D we used the self-duaHty relation ( |135| ). 

We are finally prepared to effect the spinorial translation of (EJ), which we rewrite as 



(142) 



We have 



kl 



V2 



(143) 



where we made use of the identity (|130|) . The 2- form will be useful in our manipulations 



V AB V -4/1' V B 

^ab ■— ^[a ^b] A' 



[a'- by 

— f\ t' fio^'rr^ 



V2 



^abcT\a 



c i AB 



(144) 
(145) 

(146) 



In (144) we used the covariant tetrad e^^ which is the inverse of E^^ given in (Q) and 



N 



(147) 



while in (145) we used the identity (130) and the definition (^. Wc may therefore reexpress 
1""^ using the identity (|l3l|): 

fa ^ l_^abc^^AB^^ (^48) 



Finally, substituting both (143) and (|14S|) into (142), we conclude that |2 



7^. 



CD 



i ABCTj CD^i]^ab 
AB 

XCDABl^ab : 



AB 



(149) 



where 



XCDAB ■— AB<^j CD^ 



I] ■ 



(150) 



Note that because of the symmetry and tracelessness of , xabc d is completely symmetric 



XABCD — X(ABCD)- 



(151) 
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All that remains is to show that xabcd is simply the Weyl conformal spinor '^abcd- We 
refer to the spinorial form of the curvature |3C[| 

TIaa'BB'CD = S^^^/S^s^/T?.^^'^^ 

= ABCD^A' B' + '^A'B'CD'^AB (152) 

Thus 

Tlab AB = ^ ABCo'^ab^^ + ^b^^ ^A'B'AB- (153) 



Comparing ( 149 ) and (153), we deduce that 

XABCD^A'B' — ABCD^A'B' + A' B' C D<^AB- (154) 

It follows that ^A'B'CD — 0; in other words the Ricci tensor vanishes, and 

XABCD = ^ABCD = CTi ABCT] CO'^ij- (155) 

The definitions of the Newman Penrose scalars are 





= *0000, 


(156) 




= *0001, 


(157) 




= 'Soon, 


(158) 


*3 


*0111, 


(159) 


^1^4 


= l-iiii, 


(160) 



so our claim - ( |38| ) has been verified. 

We conclude by giving the general tensorial form of the Bel-Robinson tensor T^yp^ |30) 

T^^p, I]/-4'E/^'E/^'l]/^'vI/^BCD*A'i3'C'i5'- (161) 
Thus in our tetrad basis the components are 

TiJKL — Sj^^ ^L*^ ^ABCO'i'A'B'C'D'- (162) 

We compute, for example, the Bel-Robinson super-energy density in this orthonormal basis for 
the planar cosmology considered in Section 4. 

^ XAA' xBB' f:CC' f:DD' .T, Tv; 
JOOOO = -0 ^ ABCD^ A'B'C'D' 

= i(a2 + a6 + fe2). (163) 
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